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In this paper we present the Uppsala Quantum Chemistry package (UQUANTCHEM), a new and 
versatile computational platform with capabilities ranging from simple Hartree-Fock calculations to 
state of the art First principles Extended Lagrangian Born Oppenheimer Molecular Dynamics (XL- 
BOMD) and diffusion quantum Monte Carlo (DMC). The UQUANTCHEM package is distributed 
under the general public license and can be directly downloaded from the code web-site [T]. Together 
with a presentation of the different capabilities of the uquantchem code and a more technical discus- 
sion on how these capabilities have been implemented, a presentation of the user-friendly aspect of 
the package on the basis of the large number of default settings will also be presented. Furthermore, 
since the code has been parallelized within the framework of the message passing interface (MPI) , 
the timing of some benchmark calculations are reported to illustrate how the code scales with the 
number of computational nodes for different levels of chemical theory. 



I. INTRODUCTION 

In order to complement the broad selection of quan- 
tum chemistry codes available with an "easy to use", 
development friendly and yet versatile computational 
framework, the Uppsala Quantum Chemistry package 
(UQUANTCHEM) [I] has recently been developed. The 
user-friendliness of the UQUANTCHEM package has 
been ascertained by a large set of default values for 
the computational parameters, in order for the unex- 
perienced user not to get overwhelmed by technical de- 
tails. Furthermore, thanks to the limited number of pre- 
installed computational libraries required prior to the in- 
stallation of the UQUANTCHEM code (only the hnear 
algebra package (LAPACK)[2 and the basic linear al- 
gebra subprograms (BLAS)|3j are required) the package 
is also very simple to install. The UQUATCHEM code 
has been written completely in Fortran90 and comes in 
three versions; A serial version, an openmp version and 
a MPI version. In the case of the serial and the openmp 
version, more or less generic make files are provided for 
the ifortran and gfortran[3j compilers. The MPI version 
of the UQUANTCHEM code comes with pre-constructed 
make files for five of the largest computer clusters in Swe- 
den, the Lindgren cluster(S], the Matter Cluster[6 , the 
Triolith cluster [7], the Abisko cluster [S] and the Kalkyl 
clusterpl. These makefiles can be used as templates to 
create makefiles for a broad selection of clusters. 

The wide range of capabilities of the UQUANTCHEM 
package is perhaps best illustrated by the different levels 
of chemical theory in which the electron correlation can 
be treated by UQUANTCHEM, ranging from Hartree- 
Fock and M0ller plesset second order perturbation the- 
ory (MP2)[Tn] to configuration interaction ^Tj, density 
functional theory (DFT) [lH [T3] and diffusion quantum 
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Monte Carlo (DQMC)[T4]. 

The UQUANTCHEM package provides a platform on 
which further development can easily be made, since the 
implementation of the different electronic structure tech- 
niques in UQUANTCHEM is to a large extent based on 
the text books of Szabo and Ostlund^l\ and Coofc|15j. 
and since the implementation is almost in one-to-one 
correspondence with these texts, i.e the code has been 
transparently written and well commented in reference 
to these texts. 



II. CAPABILITIES 

The UQUANTCHEM code is a versatile computa- 
tional package with a number of features useful to any 
computational chemist. The main ingredient in any 
quantum chemical calculation is the level of theory in 
which the correlation of the electrons are treated, here 
the UQUANTCHEM package is no exception. The least 
computational demanding level of theory explored by the 
UQUANTCHEM code is the Hartree-Fock level of the- 
ory, where the electron correlation is completely ignored. 

In the context of Hartree-Fock total energy calcula- 
tions it is also possible to calculate analytical inter- 
atomic forces, enabling the user to either relax the molec- 
ular structure with respect to the Hartree-Fock total 
energy, or perform molecular dynamics (MD) calcula- 
tions. Here the user can either choose to run a Born- 
Oppenheimer Molecular Dynamics calculation (BOMD), 
or an extended Lagrangian Molecular Dynamics calcula- 
tion (XL-BOMD), where the density matrix of the next 
time-step is propagated from the previous time step by 
means of an auxiliary recursion relation [I6]. The ad- 
vantage of the XL-BOMD methodology over the BOMD 
approach is that in the case of XL-BOMD, there is no 
need of a thermostat and an accompanying rescaling 
of the nuclear velocities in order to suppress any en- 
ergy drift. In Figure [T] the results obtained with the 
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FIG. 1: Total energy of a Born-Oppenheimer molecular dy- 
namics (BOMD) calculation of a H2O molecule, where the 
density matrix from the previous time step, is used as the ini- 
tial guess to the SCF optimization with the energy converged 
to < 0.01 ^Hartree, blue dotted line. The total energy of a 
Extended Lagrangian Born-Oppenheimer molecular dynam- 
ics (XL-BOMD) [16] calculation of a H2O molecule with 5 scf 
iterations per time step, full black line. 



UQUANTCHEM code for a. H2O molecule using tlie 
XL-BOMD and BOMD sctiemes wittiout tliermostat are 
shown. Here the superiority of the XL-BOMD approach 
over the BOMD scheme is manifested by the lack of en- 
ergy drift in the former's total energy. 

On the intermediate level of electron correlation theory 
implemented in the package one finds MP2[10 , CISD[TT] 
and DFT[12l|l3]. When using the DFT level of electron 
correlation it is also possible to calculate analytical inter- 
atomic forces, and therefore also perform structural re- 
laxation and molecular dynamics calculations. Here the 
DFT forces are analytical to the extent that the gradients 
with respect to nuclear coordinates of the exchange cor- 
relation energy are calculated as analytical gradients of 
the quadrature expression used to calculate the exchange 
correlation energy [17]. 

The highest level of electron correlation theory pos- 
sible to utihze within the UQUANTCHEM package is 
DQMC[T3]- Here it is possible, within the fixed node 
approximation [15] . to calculate total ground state ener- 
gies of medium sized molecules taking into account > 
90% of the correlation energy. In Figure [2] the estimated 
charge density of a H2O molecule is shown, here calcu- 
lated with DQMC as implemented in UQUANTCHEM. 



III. IMPLEMENTATION DETAILS 

The code utilizes a localized atomic basis set, where 
each basis function, 0^, is constructed from a contraction 



FIG. 2: The estimated charge density in water using a diffu- 
sion quantum Monte Carlo algorithm as implemented in the 
UQUANTCHEM code. Here red colors correspond to high 
values of the charge density and blue or bluish colors corre- 
spond to low values. The charge density was estimated from 
2000 walkers and 4.0 • 10® time steps with a time step of At = 
1.0 ■ 10~* an, and a cusp corrected[19j cc-pVTZ basis. The 
resulting ground state energy was Eo = —76.4409 ± 0.0002 
an. The rendering of the charge density map was obtained 
by using the UQUANTCHEM output file CHARGEDENS . dat as 
input to the Matlab script chargdens2DIM.m, which is a script 
provided with the UQUANTCHEM package. 



of primitive cubic gaussian orbitals[T5] 



(j)i{x,y,z) = 

n 



Here, din are the contraction coefficients, = 
{Xi, Yi, Zi) are the atomic coordinates at which the basis 
function is centered, ai are the primitive gaussian expo- 
nents and Li, Mi and Ni integer numbers determining 
the angular momentum, I = Li + Mi + Ni, of the corre- 
sponding basis function. In what follows we will suppress 
the spin part of the basis functions and assume that the 
spin degrees of freedom are treated implicitly, i.e have 
been integrated out. 

Thanks to the use of primitive gaussians almost all the 
matrices involved in the different implementations, such 
as the overlap matrix, S, the kinetic energy matrix, T 



3 



and the nuclear attraction matrix, V defined by 



Sij — 
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have been calculated analytically. Here, Z„ denotes the 
atomic numbers of the nuclei. The implementation of the 
analytic evaluation of the above integrals follows almost 
exactly the outline given in D. B. Cook's book Handbook 
of Computational Chemistry. 

In order to enhance the performance of the code, the 
electron-electron integrals. 



have been calculated by Rys quadrature [20j, even though 
they can be calculated analytically as is described in 
Cook's bookdi]. 

The exchange correlation energy, and the corre- 
sponding exchange correlation matrix elements, V^^"^ 
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defined through the exchange correlation energy density, 
£xc[p]i E^nd the functional derivative of the exchange cor- 
relation energy, have been calculated by decom- 
posing the above spatial integrals into sums of integrals 
over atom-centered "fuzzy" polyhedra, as described by 
Becke[2T]. Each one of these polyhedral integrals is com- 
puted by the use of Gauss-Chebyshev quadrature of sec- 
ond orderill], and Lebedev quadrature [23 . Here, p de- 
notes the charge density. 

The only difference between the implementation of the 
above integrals in UQUANTCHEM and the outline given 
by BeckegH, is that in the UQUANTCHEM implemen- 
tation, the mapping of the interval x £ [—1,1] into the 
radial integration interval r S [0, cxd], enabling the use 
of Gauss-Chebyshev quadrature, is contrary to what is 
prescribed by Becke done with the mapping 



r = -r,Jog 



1 - 
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Here r,„ is the Slater atomic covalent radius ^24] of the 
atom at which the corresponding "fuzzy" tetrahedron is 
centered. It has been shown by Mura et aZ[25] that the 
above mapping results in a much more accurate numeri- 
cal integration as compared to what is achieved with the 
mapping proposed by Becke |21). 



The density functionals provided with the 
UQUANTCHEM package are the local density ap- 
proximation (LDA) functional of Vosko et a/|53], the 
revised gradient corrected functional of Perdew et aZ[27] 
(revPBE)[28, and the B3LYP hybrid functional [26, i2i- 

The MP2 implementation and the CISD implementa- 
tion very closely follow the outline given in the text-book 
of Szabo and OstlundfTT]. 

The implementation of the DQMC algorithm in 
UQUANTCHEM follows closely the algorithm outlined 
in the work of Umrigar et al |14j . However, in 
UQUANTCHEM the trial function is constructed with 
a much simpler Jastrow factor, ^J, and the slater de- 
terminants are constructed from cusp corrected gaussian 
orbitals, instead from Slater type orbitals (STO), as in 
the work of Umrigar et al. The implementation of the 
cusp correction in UQUANTCHEM follows the prescrip- 
tion given by S. Manten and A. Liichow [19 . The explicit 
form of the trial function used for the importance sam- 
pling in the DQMC of UQUANTCHEM is given by: 



Where 
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fij = ~ ''il is the distance between electron i and 
J, and are the slater determinants created from 
spin up respectively spin down orbitals. The orbitals are 
constructed from the unrestricted Hartree-Fock (URHF) 
self consistent solution. Here 5 = 0.25 if the spin of the 
electrons i and j are identical otherwise if the spins are 
opposite, 5 = 0.5. Here, b and c are Jastrow parameters 
that can be chosen and optimized by the user. 



IV. TECHNICAL DETAILS 

As was mentioned in the introduction, the 
UQUANTCHEM package has been written with 
the aim of keeping a high level of transparency in order 
to facilitate further development, with the result of a 
somewhat limited computational speed for the serial 
implementation of the code. In this section we will 
mainly focus on how this limitation in computational 
speed has been dealt with by means of parallelization 
within the context of the message passing interface 
(MPI). 

In Figure [3] the normalized computational times 
for a couple of benchmark calculations of the 
UQUANTCHEM MPI version are displayed. Here the 
lowest gain in computational speed can be found for 
the Hartree-Fock type of calculations (HF), where the 
speedup is saturated at ^ 500 processors resulting only 
in a ~ 60% decrease of the computational time. Whereas 
the best speedup can be found for the MP2 and CISD 
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FIG. 3: Normalized computational time as a function of the 
number of processors used, for calculations performed at dif- 
ferent levels of theory. Here the black circles are the com- 
putational times of a Hartree-Fock calculation of H2O using 
a cc-pVQZ basis set, the red squares are the computational 
times of a density functional theory (DFT) [121 [13] calcula- 
tion of H2O using the Perdew, Burke and Ernzerhof (PBE) 
functional 27 , the green triangles are the computational times 
of a second order MoUer Plesset (MP2) ^Oj calculation of H2O 
using a cc-pVTZ basis set and the empty blue diamonds are 
the computational times of a singles and doubles configura- 
tion interaction (CISD) calculation of H2O using a cc-pVTZ 
basis set. Here to =13 s for the Hartree-Fock calculation, to 
= 68 s, for the DFT calculation, to = 592 s for the MP2 
calculation and to — 425 s for the CISD calculation. 



calculations with a ^ 95% decrease of the computational 
time, and and no saturation of the speedup can be found 
at ^^500. The difference in enhanced performance be- 
tween HF calculations and the MP2 and CISD calcula- 
tions, comes from the fact that in the case of the HF 
calculation, only the computation of two-electron inte- 
grals, {ij\kl), and the contraction of these integrals with 
the density matrices Pij 
Kij, matrices , i.e 



into Fock, J^-, and exchange. 



kl jk 

are parallelized. Whereas in the case of the MP2 calcula- 
tion, also the sum of two-particle excitations [llj are made 
in parallel, and in the case of the CISD calculations, both 
the construction of the Hamiltonian and the diagonaliza- 
tion are made in parallel, thanks to the utilization of 
the Scalable Linear Algebra Package (SCALAPCK)[32j 
divide and conquer diagonalization routine PDSYEVD 33] ■ 
Furthermore, since the storage of the CISD Hamilto- 
nian is shared amongst all the computational nodes tak- 
ing part in the calculation, memory bottlenecks can be 
avoided, permitting the computation and diagonalization 
of Hamiltonians of substantial size. 

Furthermore, in order to accelerate the convergence of 
the HF and DFT self consistent calculations, the direct 



inversion of the iterative subspace algorithm (DIIS) of 
P. Pulay[34j [35] has also been implemented in both the 
serial and parallel versions of UQUANTCHEM. 



V. USING THE CODE 



The UQANTCHEM package has been developed to be 
easily installed and run on UNIX-type platforms. The 
code has been compiled without any problems on ma- 
chines running under Mac OSX 10.7, Red-Hat Linux and 
Ubuntu, using cither the gfortran or ifort compilers. 

The interface of the code with the user goes mainly 
through one single file, the INPUTFILE-file, through 
which the user specifies what type of calculation is going 
to be performed on what type of molecule. Below, the 
INPUTFILE-file specifying a DFT calculation of the total 
energy and HOMO-LUMO orbitals of a water molecule 
is given 

CORRLEVEL B3LYP 
TOL l.OE-8 
WHDMOLUMO .T. 
Ne 10 
NATOMS 3 
ATOM 1 
ATOM 8 
ATOM 1 



0.4535 1.7512 0.0000 
0.0000 0.0000 0.0000 
-1.8090 0.0000 0.0000 



Apart from the INPUTFILE-file the user is only required 
to provide the code with the so called BASISFILE-file, 
specifying which type of gaussian basis set to use in 
the calculation. The most common basis-sets used in 
modern quantum chemistry are provided within the 
UQUANTCHEM package. There is also the possibility 
to use the information at the Basis Exchange Portal, 
located at https : //bse . pnl . gov/bse/portal] to 



extend the default library of basis-sets that comes with 
the UQUANTCHEM package. 

Once the user has made sure that the files INPUTFILE 
and BASISFILE are located in the same directory as the 
code is being executed, then, in order to run for instance 
the serial version of the code, it is just a matter of giving 
the command: 

" commad-prompt> . /uquantchem . s" 
on the command line. 

Thus from the above given expose, it becomes appar- 
ent that the user can, after the code has been compiled, 
quickly start performing quantum chemical calculations 
without getting overwhelmed by a plethora of input- 
parameters and input-files. However, when necessary, 
the user can always bypass the default settings and more 
carefully specify the calculation with the computational 
input-parameters that were initially "hidden". 
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FIG. 4: The highest occupied molecular orbital (HOMO) 
iso-surface, (a), and the lowest unoccupied molecular orbital 
(LUMO) iso-surface, (b), of the Alanine amino acid calculated 
with UQUANTCHEM using a 6-31G** basis set together 
with the B3LYP [261I29H3T] functional. Here the rendering of 
the iso-surfaces was obtained by using the UQUANTCHEM 
output files HOMO . xsf and LUMO . xsf as input to the xcrysden 
graphical software package |36| . 



VI. INTERFACING WITH OTHER 
COMPUTATIONAL SOFTWARE 



The UQUANTCHEM software comes with a set of 
supporting utility scripts for generating two-dimensional 
charge-density plots and three-dimensional iso-density 
plots with Matlab. See for example Figure [2] showing a 
two-dimensional charge-density plot of a water molecule 
generated with the chargdens2DIM.m Matlab script in- 
cluded in the software package. 

The UQUANTCHEM has been adapted to work with 
the graphical software package xcrysden |36j . by gener- 
ating a number of files in the xsf- format. In Figure |4] the 
highest occupied molecular orbital (HOMO) iso-surface 
and the lowest unoccupied molecular orbital (LUMO) 
iso-surface of Alanine are shown. The rendering of the 
iso-surfaces was obtained by using the UQUANTCHEM 
output files HOMO. xsf and LUMO. xsf as input to the 
xcrysden [36, graphical software package. 
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